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The CoGeNT experiment, dedicated to direct detection of dark matter, has recently released excess events that could be interpreted 
as elastic collisions of -10 GeV dark matter particles, which might simultaneously explain the still mysterious DAMA/LIBRA mod- 
ulation signals, while in conflict with results from other experiments such as CDMS, XENON- 100 and SIMPLE. It was shown that 
5-15 GeV singlino-like dark matter candidates arising in singlet extensions of minimal supersymmetric scenarios can fit these data; 
annihilation then mostly proceeds into light singlet-dominated Higgs (pseudo)scalar fields. We develop an eff'ective Lagrangian 
approach to confront these models with the existing data on cosmic-ray antiprotons, including the latest PAMELA data. Focusing 
on a parameter space consistent with the CoGeNT region, we show that the predicted antiproton flux is generically in tension with 
the data whenever the produced (pseudo)scalars can decay into quarks energetic enough to produce antiprotons, provided the anni- 
hilation S-wave is significant at freeze out in the early universe. In this regime, a bound on the singlino annihilation cross section 
is obtained, (crv) < 10"^^ cm^/s, assuming a dynamically constrained halo density profile with a local value of p© = 0.4 GeV/cm^. 
Finally, we provide indications on how PAMELA or AMS-02 could further constrain or detect those configurations producing 
antiprotons which are not yet excluded. 
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1. Introduction 

The CoGeNT 111 Collaboration has recently claimed for ex- 
cess events in their data, while not confirmed by complemen- 
tary searches like XENON-100 |2,^ or SIMPLE L4J. Although 
CDMS-II had reported a few excess events O, while not statis- 
tically significant, their latest results are also in conflict with the 
CoGeNT data [6J . The situation has become even more confus- 
ing, though quite exciting, since CoGeNT announced hints for 
annual modulation of the event counting rate |7|. Scatterings 
of weakly interacting massive particles (WIMPs) with masses 
around 10 GeV ofl' the detector targets could explain these mea- 
surements, and interpretations in some specific particle physics 
models have been discussed in e.g. [Tsll and |9| (see also |10| 
for less conventional models). Intriguingly, this mass range 
turns out to be close to what is needed to fit the still mysteri- 
ous DAMA/LIBRA signals lEJIElini. Nevertheless, the au- 
thors of Ref . 1^ have notably shown that some of the parameter 
space individually favored by the cited experiments may not be 
compatible among each other (see also |[T4|[T5l[T6ll ). 



Email addresses: davidg. cerdenoOuam. es (David G. Cerdeno), 
timur . delahayeOuam . es (Timur Delahaye), Iavalle@in2p3 . f r (Julien 
Lavalle) 

^Multidark fellow 



These controversial data have triggered a lot of activity in the 
field, and diff'erent phenomenological directions have been pur- 
sued. While light neutralinos arising from the minimal super- 
symmetric extension of the standard model (MSSM) are still 
debated (see arguments against this proposal in e.g. ifTTl . to 
which very detailed highlights are opposed in (Wl), singlet 
extensions of supersymmetric scenarios (including the next- 
to-MSSM — NMSSM) have received a particular attention, 
e.g. |[l9l[20l|2l]|22l[2l. The NMSSM, which provides a very 
elegant way to solve the so-called yu-problem occurring in the 
MSSM (see ||24l[25l for reviews), could actually lead to a phe- 
nomenology consistent with the direct detection hints, though 
not generically, and some interesting regions of the parameter 
space are worth being investigated more deeply 1261 [211 1221 . 
In slight contrast, more general singlet extensions uncorrected 
with the yu-problem can quite easily be tuned to accommodate 
the data, e.g. |[23]| , because more flexible. In these models, the 
dark matter candidate is usually a light singlino-dominated neu- 
tralino. 

Interesting constraints on light WIMPs may actually come 
from colliders | 27 , 28 ,29], but most of the singlet models men- 
tioned above seem to escape them ||2T]|22l. From the astrophys- 
ical point of view, the annihilation of such light WIMPs may 
also generate gamma-ray fluxes at the edge of being excluded 
with current measurements 1301 ; observational limits on high- 
energy neutrinos from the Sun also give interesting constraints. 
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as shown in e.g. [[3T|. Another very powerful astrophysical con- 
straint comes from measurements of local cosmic-ray antipro- 
tons, which severely limits direct annihilation into quarks |32|. 
In this paper, we aim at extending the latter analysis to singlet- 
like phenomenologies in which singlino-like dark matter parti- 
cles;^ annihilate into light scalar and/or pseudo-scalar particles 
(generically referred to h and a in the following), which can 
further decay, for instance into quarks. For the sake of gener- 
ality, we will adopt a model-independent approach by using an 
effective Lagrangian relevant to any particular singlet model, 
and will derive constraints on all involved masses and cou- 
plings. We note that indirect detection constraints come into 
play only when annihilation is not velocity- suppressed, which 
implies that the dominant annihilation channel is XX ~^ ^h. 
This defines the mass range we are going to focus on, namely 
2m^ > rria -\- Jfih. 

The outline of the paper is the following. In Sect. |2] we 
present our effective Lagrangian approach and discuss the an- 
nihilation cross section and issues related to the relic den- 
sity. In Sect. [3] we briefly review the calculation of the spin- 
independent WIMP-nucleus cross section, and delineate the pa- 
rameter space to be considered for interpretations of the direct 
detection hints. In Sect.|4j we derive the cosmic-ray antiproton 
spectra generated by decays of light scalar and pseudo-scalar 
Higgs bosons, and we determine the constraints current flux 
measurements put on the effective parameter space. We con- 
clude in Sect.|5] Further details on our calculations are given in 
the appendices. 



2. Effective Lagrangian approach to singlino-like dark 
matter scenarios 



The phenomenology of light singlino-like dark matter sce- 
narios is mostly set by (i) the couplings between the neutralino, 
a Majorana fermion, and the light scalar and pseudo- scalar par- 
ticles of the Higgs sector, (ii) the masses of the involved fields, 
and (iii) the couplings of the light scalars to standard model 
particles. We address the former point here; the latter will be 



discussed in Sect. 2.2 Denoting;^ the neutralino field and 0/ a 
generic field that contains a scalar part hi and a pseudo-scalar 
part Qi such that (pi = hi -h iai, we consider the following effec- 
tive interaction Lagrangian: 



i ij,k 



Mjk_ 

nijk 



pi(Pj(Pk 



-h h.c. 



(1) 



where we Aij^k is a dimensional real coupling (rjijk is a combi- 
natory factor), and where 



Cy 



= c^i-^-ic^iYs =C^ 



(2) 



features both the couplings of the scalar and pseudo-scalar 
fields (pi to the dark matter fermionic field x- We further de- 
fine c^i (scalar coupling) and c^i (pseudo- scalar coupling) as 
real numbers. 



2.1. Higgs sector 

We consider singlet extensions of the MSSM in which the 
lightest neutralino is singlino dominated, and the light Higgs 
scalar and pseudo-scalar are singlet dominated. Nevertheless, 
the physical Higgs fields are still mixtures of the Higgs doublets 
Hu and H^ and of the singlet field S , which allows us to define 
the mixing matrix S as 



Se s^ 

—Cq S(p — Sq C^i^ C^p 
Cq C^^ - Sq C<^ 



-Cq S^ 
Cq Cif^ Ctp — Sq Siff 
Cq Si^ C(p — Cq Ci^ 



Cij/ S(p 
S\f/ S(p 



(3) 



where Sk and Ck stand for sin^ and cos^, respectively. The 
mass eigenstates corresponding to the lightest pseudo-scalar 
and scalar fields, which we respectively denote a and /z, are 
conventionally given by the first row: 



h 



s^ [sq Rc(Hd) - Cq Rc(Hu)] + c^ Rc(S) 



(4) 



a = s^^ [sQa lm(Hd) - CQa Im(if«)] + c^^ \m{S ) 



= s^^ [Amssm] + c^^ lm(S) , 



where Amssm corresponds to the same linear combination of the 
imaginary parts of Hu and H^ that describes the pseudo-scalar 
Higgs in the MSSM. The mixing angles can be determined in 
specific models from the mass terms and couplings featuring 
the superpotential and the soft potential. The Higgs mixing will 
actually play an important role not only when calculating the 
direct detection signals, as discussed in Sect. [3] but also to de- 
termine the decay channels of the light Higgs fields, which the 
antiproton production depends on (see Sect. |4.1[ ). It also helps 
us understanding the collider constraints, as discussed later on. 

In this paper, we are assuming the lightest CP-even Higgs 
to be singlet-like, in order to avoid collider constraints, and 
therefore we take c<^ > 0.9. For the two heavy CP-even Higgs 
bosons, we adopt a MSSM-like scenario in which the lightest 
one is of w-type (Hu) and the heavy one is approximately //j. 
Thus, we fix Cq = Ci^ = 0.9. 

Last but not least, as it is customary in the MSSM, we define 
tan JS as the ratio of the Higgs doublet fields vacuum expectation 
values, tanyS = {Hu)l{Hd). 

2.2. Couplings to standard model fermions 

To complement this effective setup, we also need to define 
the couplings between the scalar and pseudo-scalar fields to the 
standard model (SM) fermions, which will further determine 
the branching ratios of their decay products (annihilation final 
states) and the strength of the WIMP-nucleon interaction (di- 
rect detection signals). For a field (pi having both a scalar and 
a pseudo-scalar part, the couplings to SM fermion fields dif- 
fers for up-type quarks (or neutrinos) and down-type quarks (or 
charged leptons) as: 



^0/,up,up 
^0/, down, down 



2Mw Sj3 

gm 



S(p Cq I 75 S(p 

cea] (5) 



S^SQ-iysS^SQa] , 



2Mwc^ 

where mu (m^) is the up-type (down-type) quark mass. 
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These coupHngs are useful to compute the decay branching 
ratios of the light Higgs fields, which will further determine 
the induced antiproton spectrum. In practice, we use the light 
(pseudo-) scalar decay widths that we derive as: 



Nc —z m^. 



2m/ 



3/2 



(6) 



4^ 



32 TT m^. 



1- 



2m0. 



where A^^ is a color factor (3 for quarks, 1 for leptons), and 
the latter expression is relevant to the cases where the process 
h ^ aa adds up to decays into SM fermions (ruh > 2 m^). 

In the following, we will only consider those interactions 
between singlinos and SM particles which are mediated by 
the scalar and pseudo-scalar Higgs bosons, assuming all other 
contributions to be negligible. In other words, we assume 
that other sectors of the complete underlying theory decouple. 
This is a reasonable assumption given the current mass lim- 
its on squarks (331 [33| and charged sleptons |35|. Moreover, 
as regards the Higgs sector, we impose that the interactions 
with standard matter are driven by light field exchanges only 
(see Sect. [3]). 

2.3. Some collider constraints 

It is important to derive ranges for our eff'ective model param- 
eters, namely masses, mixing angles and couplings, which are 
consistent with current collider constraints. These constraints 
necessarily involve the couplings to SM particles, which have 
been made explicit in Eq. Q. We basically use the constraints 
provided in |23 | for light scalars, and in |36| for light pseudo- 
scalars. The only additional feature we impose to our singlino 
dark matter candidate is to reproduce the region compatible 
with the WIMP interpretation of the CoGeNT signal, given 
the parameter space delineated for the scalar and pseudo-scalar 
fields below. 

Limits can be derived because of the potential contributions 
light scalars and pseudo-scalars may give to rare mesons de- 
cays. From measurements of B decays by Belle |37 | and SM 
background calculations by [38], the following limit has been 
derived by the authors of ll23ll for the light scalar mass, given its 
non- singlet content s^^: 



Br(B ^h-\-Xs)x Br(/z ^ yu^ yu") = 



(7) 



0.058 



0.1 



ml 



X Br(/z ^ yu+ -h yu") < 2.5 X 10" 



The decay branching ratio Br(/z ^ yu^ -hyu") is readily computed 
from Eq. ([5]). This roughly sets a lower limit for nih in terms of 
the bottom quark mass, namely nth ^rnt. 

There is an additional and independent limit on s^p which can 
be derived from LEP data on Higgs searches, especially from 
the process e'^ e~ Zh vvh \39\. This yields the upper 
limit s^p < 0.08, valid as long as mh < 2ma, i.e. as long as the 
decay channel h ^ aa is forbidden [231 . 



Constraints on the singlet content of the pseudo-scalar can 
be extracted from e.g. the CLEO Higgs search data from the 
T(15') decays |40|. They are expressed in terms of the reduced 
coupling Xd of the light pseudo- scalar a to down-type quarks 
and leptons, which is normalized with respect to the coupling 
of the CP-even SM Higgs boson, and a very conservative limit 
is given by: 



Xd = s^^ tanyS < 0.5 



(8) 



as read off' from Fig. 1 of Ref. f36l. This is definitely not a strin- 
gent limit for the phenomenology under scrutiny here, since this 
has no consequence on the pseudo-scalar decay products (the 
mixing angles factorize out in the decay branching ratios). This 
still ensures that the singlino annihilation into SM fermions via 
light pseudo-scalar exchange is always suppressed with respect 
to the annihilation into light (pseudo-)scalars. 

Finally, we set a lower limit of m^ > 1 GeV in the following 
to avoid relying on too large fine-tuning zones in more theoret- 
ically constrained parameter spaces. 

2.4. Relic density 

One of the most important constraints on dark matter mod- 
els obviously comes from requiring that they lead to the cor- 
rect cosmological abundance, which usually severely reduces 
the available parameter space. Despite the uncertainties af- 
fecting our knowledge of the expansion rate of the universe 
before the big-bang nucleosynthesis (BBN) (see e.g. f4T\), 
the general formalism to derive the relic abundance, based on 
solving a covariant Boltzmann equation, has been established 
and tested for a long time |42 lEl ISl |45| ggl gll. There are 
some public codes like DarkSUSY EHl or MicrOMEGAs B9ll , 
which numerically solve the Boltzmann equation, originally de- 
signed for supersymmetric models. Both these codes essen- 
tially make use of the same numerical scheme and of the same 
tables for the effective degrees of freedom, which rely on earlier 
works B4l |46l |47]| . Consequently, diff'erences between them 
are expected to come from the annihilation cross section calcu- 
lation only. An alternative to these numerical tools, but at the 
cost of loose accuracy, is to use approximate analytical expres- 
sions ||45l|46l|50l. Here, we use our own numerical tools, more 



detailed in Appendix B to which we refer the reader for the 
definitions of the relevant parameters. Denoting 7>/ the comov- 
ing dark matter density after decoupling such that 7>/ » Yq^, 
its value today Yq can be calculated from the following expres- 
sion: 



Yn = 




(9) 



where x = m^/T (T is the thermal plasma temperature). Note 
that the second term in the right-hand side, featuring the an- 
nihilation cross section (crv) (see Sect. |2.5| ) and the eff'ective 
degrees of freedom g^, is usually dominant. A numerical code 
is basically useful to compute y>/ = F(x>/) to a good accuracy 
from the full Boltzmann equation (see Eq. |B.1| in [Appendix] 
[BJ, given the temperature-dependent annihilation cross section 




Figure 1 : Left: effective degrees of freedom relevant to the relic density calculation, where a sharp first order QCD phase transition is assumed. Right: impact of 
assuming different QCD phase transition temperature on the relic abundance of Majorana fermions — we fix the annihilation cross-section to 3 x 10~^^cm^/s for 
the sake of illustration. 



(crv) and the effective degrees of freedom g^. The reUc abun- 
dance then merely reads: 



soYoh' 



^ Pc ^ 3x1002 

^ 2.144 X 10^ Yo , (10) 

where we have translated the critical density pc in terms of the 
Hubble parameter H, using the relation = (100 x /z)^ = 
87rpc/3. The parameter sq ^ 2891.23 cm"^ is the entropy den- 
sity today — we have taken a CMB temperature of 2.72548° K 
EH. 

One of the key issues in calculating the relic density of 
-10 GeV WIMPs is the QCD phase transition. Indeed, these 
WIMPs typically decouple at temperatures Tf ~ ~ 400 

MeV, when quarks (and gluons) go from asymptotic freedom 
to confinement into hadrons. This transition is therefore char- 
acterized by an important change in the relativistic degrees 
of freedom appearing in Eq. ([9]), which can be a source 
of systematic errors in the predictions (as was already clear 
from e.g. Figs. 5 and 6 of |44|). For instance, slightly mod- 
ifying the phase transition temperature Tqqd can lead to sig- 
nificant changes in the final relic abundance, by up to 50% (a 
factor of ~ 4 in g^), for those models which happen to decou- 
ple after instead of before the transition, after modification of 
Tqcd (see Eq.|9]). The point is that we do not know much about 
this phase transition yet despite the recent progresses obtained 
in lattice QCD |52|, and theoretical uncertainties in this regime 
might still be sizable 1531 . The current state of the art |[54|| is im- 
plemented in e.g. DarkSUSY, but modifications of the transition 
temperature is not user-free. 

To try to bracket potential theoretical uncertainties coming 
from this phase transition, we have recomputed the eff'ective 



degrees of freedom for diff'erent values of Tqcd, considering 
a close to first order phase transition. An example is shown in 
the left panel Fig.jl] where we display geff^^^, ^eff^^^ and g^^ 
for two diff'erent assumptions Tqcd = 150 or 400 MeV (see 
Eq. B.2 in [Appendix B for the definitions). The right panel 
exhibits the consequence on the dark matter relic abundance, 
assuming a constant annihilation cross section set to its canoni- 
cal value (crv) = 3x 10"^^ cm^s"^ and only varying the WIMP 
mass m^. We clearly see that in the region of interest, - 10 
GeV, the relic density can decrease or increase by a factor of 
~ 2, depending on whether decoupling occurs just before or 
just after the QCD phase transition — this can be understood 
from Eq. ([9]), where we see that Fq ^ (sH^ (<^v))"^ Therefore, 
some models in this mass range can easily be cosmologically 
excluded for one choice of Tqqd but allowed for other choices. 
In the plot, as will be the case in the rest of the paper, we have 
used the WMAP-7 data for the present dark matter abundance, 
i.e. Vl^h^ = 0.1109 ± 0.0056 [55], which appears as the green 
band. 

This illustrates why uncertainties on the quark-hadron tran- 
sition and on Tqcd should be taken into account, as often un- 
derlined in the past, like in 1561 . A detailed implementation of 
this transition is beyond the scope of this paper. Still, we con- 
sider the two extreme cases depicted in the left panel of Fig.[T] 
which is likely sufficient to encompass the full theoretical error. 
Finally, we stress that the impact on indirect detection signals 
amounts to the same factor of ~ 2. The annihilation cross sec- 
tion needs to be relatively larger to get the correct abundance 
when dark matter particles decouple after the QCD phase tran- 
sition, as clearly seen in the right panel of Fig. [T] We have 
checked that using other predictions for the evolution of the rel- 
ativistic degrees of freedom and their change during the QCD 
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phase transition, for instance those included in the DarkSUSY 
package, led to a < 10% difference in the relic density predic- 
tions, or equivalently in the annihilation cross sectioij^ There- 
fore, taking a 3-cr band for ^^h^ is well enough to encompass 
the theoretical uncertainties in the effective degrees of freedom 
for a given Tqcd- Thus, we will consider any configuration as 
cosmologically relevant whenever leading to a relic density of 
Q;^/z2 = 0.1109 + 0.0168. 

2.5. Annihilation cross section 

The annihilation cross section is relevant not only to indirect 
detection, but also to direct detection, since in both cases it de- 
termines the cosmological abundance today. In the following, 
we assume that the only particles involved in the annihilation 
process are the light scalar and pseudo-scalar Higgs bosons, h 
and a respectively. This is justified by the fact that we are in 
the very light mass and singlet/singlino dominated regime — 
other Higgs bosons, charginos and sfermions have much larger 
masses, while similar couplings. In such a regime, annihilation 
implies the following final states: XX cici, XX hh and 
XX cih, which can be summarized thanks to the Feynman 
diagrams represented in Fig. [2] 

As it is well known (e.g. 153), xX cih is the only tree- 
level process among the three above possibilities to have a non- 
vanishing zero- velocity limit, i.e. the only process relevant to 
indirect detection. It is kinematically allowed as long as 2 m^^^ > 
nia + Mh, which defines the mass range we are dealing with 
in this analysis. The full calculation of the annihilation cross 

Here, we just summarize a 



section is detailed in Appendix A 



few helpful results. It usually proves interesting to expand the 
thermally averaged annihilation cross section as a power series 
of the inverse of the variable x = m^/T oc (xf ~ 20 at freeze 
out, and ~ 10^ in the Galaxy today): 

{crv) = a-^b/x-^0(x-^) . (11) 

Both coefficients a and b take analytical expressions in the 
present study. The S-wave contribution is given by (we neglect 
the widths of the exchanged particle^: 



(1- 



2 m 



'X 



)(1- 



2 My 



) 



(12) 



This equation fully sets the annihilation cross section for the 
indirect detection signals that we will consider later. We im- 
mediately see that there is a competition between the s and t/u 
diagrams through the couplings Aaah and c^^h- Note that the lat- 
ter is also the one that fixes the contributions of CP-even Higgs 
bosons to the direct detection signal (more details in Sect. [3]), 
which sketches a way to link the indirect detection constraints 
to the direct detection signal: the relation is straightforward in 
the case Aaah = 0, corresponding not only to the suppression of 
the ^-channel of ah but also of that of process aa. 

Finally, for Aaah ~ c^xh close-to-degenerate masses, the t/u 
diagram dominates the cross section, which again induces a sig- 
nificant correlation between the direct detection signal and the 
indirect detection constraints. 

We also remark that the coupling c^^a does drive the a-term 
amplitude, which is actually quite obvious since xX is bound 
to be a pseudo-scalar state at null velocity. Suppressing this 
coupling is a very simple way to avoid indirect detection con- 
straints. In that case, annihilation before freeze-out would 
proceed mostly through xX hh, but also partly through 
XX (^CL via the ^-channel, provided these final states are kine- 
matically allowed. 

Although the connection between the elastic interaction with 
nuclei and the annihilation into ah may be understood from 
the above equation, we still have to make sure that the a-term 
of the cross section expansion is the dominant one in the early 
universe, in order for it to be fixed by the relic abundance. In- 
deed, should the processes xX ~^ hh be dominant at freeze 
out, then that would lessen the bounding strength of indirect 
detection and also make us lose potential correlations between 
direct and indirect signals. We recall that the correlation is def- 
initely lost when > ruh (ma) and 2m^ < ma+ mu, for which 
XX hh (aa) is the only annihilation possibility left; we do not 
consider such cases here. 

For completeness, we provide the first order terms baa and 
bhh associated with annihilation into a a and h h, which drive the 
relic density when 2m^ <ma+mh (still neglecting the particle 
widths): 



XXa 



4m^ 



^aah my 1 



(ml + ml) 
4m^ 



^We have also made some comparisons between our calculation and the one 
provided in MicrOMEGAs, and found an agreement within ~ 1% when using 
the MicrOMEGAS table for the degrees of freedom (identical to the DarkSUSY 
table). 

^Particle widths are included in the full numerical calculation. 



4 



1 - ml/m^ 



(13) 



16c' <K-<r + 



(mr — m ) 

8 Aaah C^xh ~4xa ^X ^x " 77^2 ^ 

(4m^-mf) 



+ 3 (Aaah Cxxh)' 



(4 ml 
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Figure 2: Feynman diagrams relevant to the^^ annih 



Ck 



^hh 



6A7i{2ml-mlf 



\6c'ml{9m'^2ml-^miml)- 



(14) 



(5m^ - 2ml) 
{Ami ~^h^ 



2\4 



3 A 



hhh 



(Am^-ml)^ 



These terms are the leading terms of the annihilation cross sec- 
tion when 2m^ < ma -\- mh. We emphasize that even in the case 
2m^ > ma -\- mh, i.e. when Gah is not suppressed, the velocity- 
dependent terms can not only partly contribute to the full an- 
nihilation cross section at decoupling in the early universe (up 
to -10-20% when couplings and masses are taken similar, and 
considering also bah, not given here), but can also dominate in 
some cases, e.g. when the coupling c^^a is suppressed. This 
means that the full annihilation cross section needs to be consid- 
ered for relic density calculations, which we have done in this 
study. In fact, we did not use this Taylor expansion approach 
to compute the thermal average of crv in the relic density cal- 
culation, but the exact numerical method presented in [46, 4T|. 
Still, the velocity independent term Gah given in Eq. ([12]) fully 
characterizes the annihilation cross section relevant to indirect 
detection. Moreover, Eqs. ([13]) and ( [T4| ) will prove very useful 
to interpret the antiproton constraints (see Sect.|4]). 

3. Direct detection of singlino-like particles 

In this paper, we aim at confronting the light singlino param- 
eter space compatible with the CoGeNT region to the antipro- 
ton constraints. It is therefore important to provide some details 
on how the CoGeNT region can be covered in the frame of our 
effective approach. 

The singlino couplings to the Higgs mass eigenstates will be 
denoted as c^j, in consistency with the Lagrangian introduced 
in Eq. ([T). 

Elastic collision of a singlino off a target nucleus would take 
place through the exchange of a CP-even Higgs along a t- 
channel, as depicted in Fig. [3] This scattering can be described 
in terms of an effective Lagrangian 157] , 



£eff D aq^xxmi 



Figure 3: Diagram describing the elastic interaction of a singlino with a quark. 

where the coupling aq. results from adding the contributions of 
the exchanges of the three CP-even Higgs bosons and reads. 



.^^ 2MwBi ml. 



(16) 



where / = 1, 2 for up(down)-type quarks, Bi = sin/S and B2 = 
cosjS. The mixing matrix S has been defined in Eq. ([5]). Note 
that the term of the above sum corresponding to our light Higgs 
boson h can be straightforwardly rewritten as a^. = c^j Chqq/^l, 
where the coupling Chqq has been defined in Eq. (|5]). 

The total spin-independent WIMP-proton scattering cross 
section yields 



SI 

x-p 



ml ml 



- f 



71 (mp -\- m^)^ ' 
where mp is the proton mass and 

fp 



qi-u,d,s 



= y fp ^ + 1 fp y ^ 



(17) 



(18) 



qi=c,b,t 



The quantities /^^ and /^^ are the hadronic matrix elements 
which parameterize the quark content of the proton (the second 
term is due to the interaction of the WIMP with the gluon scalar 
density in the nucleon). We will use = 0.023, /J = 0.033 
and fs =0.26 which follows from the latest results on the pion- 
nucleon sigma-term [ 58 1, which affects the determination of the 
strange content of the quark. 

The leading term in the scattering cross section is usually due 
to the coupling to the s quark (the contribution from quarks u 
and d being almost negligible and that from /^^ being at least 
a factor 5 smaller in the determination of P). Under this as- 



(15) sumption it is possible to derive from Eq. ( 16) a condition for 



the lightest Higgs contribution to be the dominant one in terms 
of the Higgs masses and the various coupHngs. In particular, 
we will demand 



and 
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(19) 



(20) 



Notice that since we are considering very small masses for the 
lightest Higgs boson, these conditions are normally easy to ful- 
fil. In our case they will imply a constraint on the Higgs doublet 
contribution and on the neutralino-Higgs coupling. For exam- 
ple, for a typical configuration with m//o = 10 GeV and mn^ = 
120 GeV the above condition implies c^i <S| > 7 x I0~^c^2 <^2- 
Under these conditions the spin-independent cross-section can 
be approximately expressed as 
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(21) 



In order to explore the region consistent with the latest hints 
for very light WIMPs in direct detection experiments [ilTimiTll. 
we have constrained the input parameters to produce a singlino- 
proton spin-independent cross section between 3 x 10"^ pb and 
3 X 10"^ pb, and singlino masses in the range 6-15 GeV. This 
window in the (cr^Lp,m^) plane contains the region compati- 
ble with the first observation by the CoGeNT collaboration [ 1 1 
and is extended in order to account for variations of this region 
when astrophysical uncertainties are taken into account O |59l 
and to include the area consistent with DAMA/LIBRA result if 
interpreted in terms of WIMP recoils on Nal (assuming negligi- 
ble channeling). Of course, the region consistent with the latest 
data release by CoGeNT 1 7 1 is also contained in the whole area. 

To try to delineate a parameter space as complete as possible, 
we have performed diff'erent random scans over the full param- 
eter space defined by the involved free masses and couplings, 
taking diff'erent values of tanyS for each. Since we require a 
given range of spin-independent cross-section and, as we said 
above, we consider only those cases where the lightest Higgs 
contribution to cr^L^ is dominant, this implies a constraint on 
the lightest Higgs composition and its coupling to singlinos. In 

leads to 



particular, for m^^^ = 6 - 15 GeV, Eq. (21 
3 X 10-4 GeV-2 < 



< 9 X 10-4 GeV 



-2 



(22) 



We note that this condition implies that for large values of tan/5 
the necessary value of the coupling c^i in order to reproduce 
the CoGeNT data can be smaller. 

In practice, the various input parameters have been randomly 
drawn in the following ranges: 

c [6, 15] GeV (23) 

c [5x 10-^5x 10-^] 

c [2C^,,,^,5x 10-1] GeV 

= 10 and 50 , 



\CxU c^i. c^?>\ 

\^ijk\ 

tanyS 



where the scan in the couplings is performed in the logarithmic 
scale. Note that the lower bound on the couplings A is not a 
fixed number, but instead twice the coupling of any light boson 
to b quarks (see Eq. [5]). This is to ensure that (the s-channel) 
annihilation dominantly proceeds into 2-body combinations of 
bosons a and h, the pair production of SM fermions being then 
suppressed and irrelevant in the computation of the relic density 
as well as of the indirect detection signals. This will induce a 
correlation between couplings X and C^h through the range im- 
posed on the spin-independent cross section [see Eqs. Q, ([16]), 
and ( [2T] )] . The lightest CP-even and CP-odd Higgs masses are 
assigned a random value, imposing the constraint mu + rua < 
2 in order to guarantee that annihilation into a pair /z, a is 
possible. We have imposed the experimental constraints on the 



Higgs masses and couplings as detailed in Sect. |2.3| and the con- 
dition that the singlino-proton spin-independent cross section 
falls within the window cr^Lp c [3 x 10-^, 3 x 10-^] pb. We have 
thus kept those configurations which fall in a region that con- 
tains the area compatible with the CoGeNT excess. As regards 
tanyS, after having tried some random scans, we have selected 
two extreme values only, which are sufficient to fully illustrate 
the impact of this parameter. In Fig.|4j we show our scan results 
in terms of as a function of for both values of tan and 
for two diff'erent QCD phase transition, with Tqcd = 150 and 
400 MeV The color index refers to relic density and antipro- 
ton constraints, and will be specified and discussed in Sect. |4. 3 1 
All samples we consider in the following are made of ~ 10^ 
random models, such as those reported in Fig.|4] 

To conclude this section, we must emphasize that some un- 
certainties in the hadronic matrix elements may aff'ect the pre- 
dictions for the direct detection rate of singlinos. The largest 
eff'ect is due to the indetermination in the strange quark content 
of nucleons. This propagates into the theoretical predictions for 
the spin-independent cross section and can be responsible for a 
variation of about an order of magnitude |[60l|6Tl|62j|63l. This 
induces an uncertainty in the values of the couplings for which 
the CoGeNT region can be reproduced and therefore enhances 
the range of these with respect to those plotted in Fig.|4] A re- 
cent illustration has been given in e.g. Fig. 5 of Ref. 11641 . We 
will further discuss this issue in Sect. 14.3.21 

4. Constraints from local cosmic-ray antiprotons 

Cosmic-ray antiprotons were proposed as a potential tracers 
of dark matter annihilation in {6^\. Dark matter-induced an- 
tiprotons have been studied in the context of supersymmetry 
|[l6||67||Ml|69l[70l|7ll and extra-dimensions Ell HD, and re- 
lated signals are likely very difficult to observe. Current obser- 
vations by e.g. the PAMELA |73 ] or BESS 11741 1751 [76l [771 
experiments are in agreement with the predictions of secondary 
antiprotons ||78l[79j|73l, i.e. those created from nuclear interac- 
tions of standard astrophysical cosmic rays with the interstellar 
gas, although some models still fail to saturate the data, leav- 
ing a small room for exotic primary contributions fSOl lSTIl . The 
case for light supersymmetric neutralinos was studied in 1821 . 

Based on recent studies on cosmic-ray transport parame- 
ters |[83l and on the local dark matter density |[84j|85l, it was 
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Figure 4: Spin-independent cross section as a function of the singlino mass, for different values of tany5 — 10 (left panel) and 50 (right panel) . Each panel shows 
-10^ models randomly drawn according to Eq. \23\ and selected such that cr^^ falls in the CoGeNT region. The color index is specified in Sect.|4.3.l] 



sliown tliat ~ 10 GeV liglit WIMPs anniliilating into quarlc pairs 
were in serious tension witli tlie low energy cosmic-ray antipro- 
ton data |32|. In tliis paper, we are focusing on a quite differ- 
ent plienomenology wliere WIMPs anniliilate into liglit scalar 
h and pseudo-scalar a Higgs bosons, which may further decay 
into quarks. These light bosons do not have necessarily the 
same masses, so one can for instance be more Lorentz boosted 
than the other. Therefore, the energy distribution of the poten- 
tially created antiprotons is expected to differ from those aris- 
ing from direct annihilation into quark pairs, and the results 
obtained in 1 32 1 can hardly be extrapolated. In this section, we 
first determine the antiproton energy spectrum associated with 
the decay of light Higgs bosons into quark pairs before calcu- 
lating the dark matter-induced flux at the Earth. 



4.1. Antiproton production 

Before estimating the antiproton flux at the Earth, one needs 
to know the antiproton spectrum before propagation, i.e. the one 
arising from the decay/hadronization of the dark matter anni- 
hilation products, here the scalar h and pseudo-scalar a. The 
only decay final states relevant to antiproton production are ac- 
tually quark pairs energetic enough to produce a pair of proton- 
antiproton. This means that antiprotons can be created each 
time a light (pseudo-)scalar boson is heavier than twice the pro- 
ton mass, roughly, which sets a production threshold to > 2 
GeV. 

The first step to get the antiproton spectrum is easy: before 
the decays of our light (pseudo-)scalars, we have a mere two- 
in-two-body annihilation so all the 4-momenta are set by kine- 
matics. Assuming that dark matter particles annihilate at rest in 
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the Galactic halo frame, the energy of resulting particle 1 is 



^1 = 



4m^ + m\-m\ 



and the norm of its momentum is 



Amy 



(24) 



(25) 



where A{a^, b^, c^) = (a^ - (\b\ + \c\)^)(a^ - (\b\ - \c\)^), which are 
enough to change from the halo frame to the rest frame of parti- 
cle 1: (y,l3)=(Ei/mi,ki/Ei). In the rest frame of particle 1, the 
quarks it decays into have energy E* = mi/ 2 and momentum 

1^*1 = ,^m^/4^ -m^. One finally gets the energy of the quarks 
in the halo frame from a mere Lorentz transform: 



Eq = 



El 

2 

El 
2 



cos(e) 



^A{4ml,mlml) I ^ 



8my 



mi 



cos(e)6. 



(26) 



Thus, the energies of the quarks and anti-quarks coming from 
the decay of particle 1 are evenly distributed between ^-S and 
^ -h 6. One finally gets the probability of having an antiproton 
of energy Ep from a quark of energy Eq f(Eq,Ep) thanks to 
the PYTHIiypackage |86|. Eventually, the antiproton spectrum 
obtained after decay of the annihilation products is: 



E 

;_ 1 o ^ ^ -T- 



i=l,2 q 



T^t,f^E„Ep) + f(E^,Ep) 
26 



dE„ 



(27) 



where the first sum is done over the annihilation products 
(l,2)=(a, h), and the second sum is over all those quark flavors 
for which 2 ^ m/ is satisfied. The branching ratios Bvi q can 
be derived from the couplings to SM fermions given in Eq. ^ 
and the decay widths given in Eq. ([6]). For the quark masses, 
we use the pole masses reported in e.g. |87|. 

It is clear that the decay into the down-type b quark- antiquark 
pairs will dominate whenever permitted. When m/ < 2 m/,, the 
dominant decay channel becomes (pi r" because of the 

more favorable dependence into tanyS with respect to the up- 
type quark c. Therefore, the antiproton production is reduced 
in the range 2mj < mi <2mh, but not completely suppressed, 
since d and s quark contributions remain significant. We have 
calculated the decay branching ratios explicitly in Fig. [5] as- 
suming c^^ = Co = 0.9 and tanyS = 10 (50) on the left (right) 
panel. Those associated with the third, second and first gen- 
eration fermions are reported with solid, dashed, dotted lines, 
respectively. These curves are relevant to h and a, but we re- 
call that we have only constrained the mixing of h so far. For 
simplicity, we take the same mixing angles for a, which can be 
considered as conservative for the antiproton production — for 



example, a slight increase of the up Higgs doublet content of a 
would reduce the production of r leptons. 

Finally, we also consider the possibility that h decays into 
2 a, a channel that is open when mh > 2 ma. This channel is 
relevant when the coupling Aaah/^h ^ Ccpqq (see Eq.[6]). Given 
the range we have considered for A (see Eq. [23]), this barely 
happens in our scan. Still, we note that if this channel is open 
and dominant, then the upper limit on the mixing parameter s^p 



can be relaxed (see Sect. |2.3| ); for simplicity, we still do not 
relax it. 

An illustration of the spectrum determination is presented 
in Fig. [6j where we have used diff'erent combinations of the 
masses m^, ma and mu and applied the Lorentz boost proce- 
dure described above — we did not include h ^ aa. We have 



considered five diff'erent cases, with m^ 



12 GeV for each. 



For this work we made use of version 6.4.24 with CDF tune A. 



We first see that each time m^,m/j > 2 GeV, the antiproton 
production is rather efficient and peaks with almost the same 
amplitude, within a factor of ~ 3, with a sharp decrease. The 
brown curve (third from top at peak energy) shows the case 
m^ = ma = mh, which means that a and h both mostly decay at 
rest into b quarks. Nevertheless, these quarks are created with 
a relatively small momentum (E^ = ma,h/'^ = 6 GeV), which 
reduces the proton and antiproton production efl&ciency. The 
violet (upper) curve shows a more favorable situation in which 
ma = 12 GeV and m/^ = 4 GeV, i.e. both a and h can decay 
into quarks with momenta large enough to fragment into an- 
tiprotons efficiently — the center of mass energy is 2m^, and 
the production of r is still not strongly dominant for the h de- 
cay because very close to threshold (see Eq.[6]and Fig. [5]). The 
pink curve (second from top) characterizes an unexpected situ- 
ation where, although h and a decays dominantly into r leptons, 
the remaining ~ 30% going into s and d quarks are very effi- 
cient at producing antiprotons thanks to their large momenta. 
The orange curve (fourth from top) shows a less favorable case 
for antiproton production in which a is not massive enough, 
and the b quarks coming from the h decay have low momenta, 
which translates into a factor of ~ 2 suppression with respect 
to the previous case. Finally, the blue (lower) curve corre- 
sponds to an even worse case, too small m« and mh = 10 GeV, 
where b quarks are inefficient at producing antiprotons in ad- 
dition to threshold eff'ects (the decay into r comes into play — 
see Fig. [5]). 

These spectral diff'erences that depend on the mass combi- 
nation demonstrate that there are non-trivial ways to suppress 
the antiproton production, e.g. producing b quarks at thresh- 
old with low momenta. Nevertheless, we will further show that 
each time a or h can decay into antiprotons efficiently, the an- 
tiproton flux is always in tension with the current data, provided 
the s-wave contribution to the annihilation cross section given 
in Eq. ([12]) is dominant at freeze out. 

4.2. Antiproton flux at the Earth 

4.2.1. Dark matter halo 

To predict the antiproton flux at the Earth, we need to specify 
the dark matter distribution in the Galaxy. Many conventional 
dark matter halo functions and related parameters have been 
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Figure 5: Branching ratios for the (pseudo-) scalar <p decay into SM fermions as a function of its mass m^. Contributions of the third, second and first generation 
fermions are reported with soHd, dashed, dotted Hnes, respectively. 
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Figure 6: Antiproton energy distributions for different mass combinations. 



widely used in the literature so far, but we aim at employing 
a distribution as constrained as possible, both theoretically and 
observationally. From the theoretical point of view, there seems 
to be a consensus in using an Navarro-Frenk- White (NFW) pro- 
file |88, 89] or an Einasto profile |90, 91 1, which are found to 
provide good fits to the dark matter distributions in cosmolog- 
ical N-body simulations (see e.g. |92|). The former is mildly 
cuspy in the center, behaving like r"\ while the latter asymp- 
totically flattens towards the center. Nevertheless, we stress that 
the very central parts of the Galaxy could difl'er from what most 



of dark matter-only N-body simulations have predicted because 
baryons dominate the gravitational potential there, which might 
modify the central behavior of the distribution |93|. From the 
observational point of view, observations of low surface bright- 
ness galaxies also seem to indicate cores instead of cusps in 
these objects |94|, showing that profiles are likely not universal 
in the centers. Still, in contrast to gamma-rays, the antiproton 
signal is well known to have almost no dependence in the dis- 
tribution in the Galactic center, provided it is less cuspy than 
r"^/^ |[69j |79l |32l : this is due to diff'usion eflfects. Much more 
important is the local dark matter density, which roughly sets 
the antiproton flux. It is therefore of stronger relevance to use a 
model dynamically constrained on kinematic data. 

Recently, complementary efl'orts were undertaken to con- 
strain the local density of dark matter p© ^2 ^|, the re- 
sults of which all point to p© ^ 0.4 GeV/cm^. Here, we will use 
the spherically symmetric NFW profile as constrained in 1951 , 
which we recall hereafter: 

1 + r©/r,''^ 



p{r) = p© 



1 + rlvs 



(28) 



the variable r denoting the galactocentric radius, and setting the 
parameters to r© = 8.28 ± 0.16 kpc, p© = 0.40 ± 0.04 GeV/cm^ 
and Vs = 18 + 4.3 kpc. We emphasize that the antiproton flux 
will be sensitive mostly to variations in p©, not in the other pa- 
rameters. 

In principle, consistency imposes to consider dark matter 
substructures |S6l[97l, which can increase the antiproton flux 
at low energy ||98j[99l. While the amplification is expected to 
be small f99\, we stress that even a factor of 2 could have im- 
portant consequences in delineating the antiproton constraints, 
as will be shown later. Here, we will not consider substructures, 
which can therefore be considered as a conservative approach. 
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L 


Vc 


Va 


kpc^/Myr 


kpc 


km/s 


km/s 


0.0112 0.70 


4 


12 


52.9 



Table 1: Transport parameters associated with the best-fit model of 11031 . 



4.2.2. Antiproton transport and flux at the Earth 

Once antiprotons are injected in the Galaxy, they diffuse on 
magnetic inhomogeneities and experience different additional 
processes. At low energy, say below a few GeV, they are sensi- 
tive to Galactic winds which translate into convection upwards 
and downwards the Galactic disk and into adiabatic energy 
losses; they are also sensitive to nuclear interactions with the 
interstellar gas, which may destroy them and/or decrease their 
energy; finally, they can be reaccelerated because of the proper 
motion of the magnetic scatterers. Higher energy antiprotons 
are roughly sensitive to spatial diffusion only. The general 
steady state equation that drives their evolution in phase- space 
reads: 



-V[{k,AE)^-V,)n] (29) 



spatial current 



current in momentum space 

+ IVV 

spallation 

where TV = N(x,p) = dn/dp = jSdn/dE is their differen- 
tial density, Kxx and Kpp are respectively the diffusion coef- 
ficients in space and momentum, Vc is the convection veloc- 
ity and = Ts(x, p) the spallation rate. This equation can be 
solved numerically QOO] or semi- analytically |101, 102, 103 1. 
In this study, we have used the semi-analytical method sketched 
in e.g. ll03ll69l[79l [98 K to which we refer the reader for more 
details. For the transport parameters, we have utilized the best- 
fit model found in 11031 and still in excellent agreement with 
more recent analyses {e.g. fST]). This transport model was 
dubbed med (for median) configuration in |69|. It is a slab 
transport model, with a cylindrical diffusion volume of radius 
7? = 20 kpc, half-height L = 4 kpc, in which the diffusion coef- 
ficient obeys K^xiE) =l3Ko(K/l GYf (^ = p/q is the cosmic- 
ray rigidity). The main parameters are reported in Tab.[T] 

The source term is characterized by the specific dark matter 
model configuration: 



Q(^,E) = Q(r,E) = 



(crv) fp(r)l ^A^dec 
nty I dE 



(30) 



where (crv) = aah is the annihilation cross section in the Galaxy 



given in Eq. ( 12 ), and dN^^cldE is the antiproton spectrum aris- 



ing from the (pseudo-)scalar decay and determined in Sect. [4^1 
Assuming a Green function Q{x,p <— Xs.Ps) solution 



lowing expression: 

(^o(^) =^^Jdps d'^s p^Xs. Ps) Q(Xs, Es) 



(31) 



In practice, we use the Bessel series method presented in |[69l 
[791 , but we have checked with a couple of configurations that 
results obtained from Green function methods 1791 Ell are 
equivalent. 

4.3. Results and discussion 

We have performed a scan over the masses and couplings 
associated with Eq. ([T]) such that we keep only those configura- 
tions falling in the so-called CoGeNT region (see Sect.[3|. For 
these configurations, we have computed the annihilation cross 
section summing over the three diagrams of Fig.[2]for all possi- 
ble final states, i.e. aa, hh and ah as soon as kinematically al- 
lowed — we recall that we have imposed 2m^>ma+mh. Then, 
we have derived the relic density assuming different QCD phase 
transition temperatures, and computed the antiproton fluxes for 
models leading to the correct relic ab unda nce, i.e. 0.0941 < 
Qv/?^ < 0.1277 in this study (see Sect.|2!4|). The flux is com- 



puted assuming the transport parameters of Tab. [T] and a solar 
modulation is subsequently accounted for in the force-field ap- 
proximation II104II1Q5L with a Fisk potential of O = 600 MV, 
consistent with the average solar activity during the PAMELA 
data taking 11731 . For the secondary antiproton flux, we use 
the predictions obtained in |78| with the same transport pa- 
rameters. It is represented in the left panel of Fig. |7] together 
with a prediction using another parameter set characterized by 
a larger diffusion halo with L = 15 kpc. Both predictions are 
consistent with the current PAMELA data |73|, while the lat- 
ter transport model would lead to a significantly larger primary 
antiproton flux because integrating over much more dark matter 
annihilation products. It turns out that large diffusion halo mod- 
els, i.e. with L > 5 kpc, are favored by several recent studies on 
Galactic cosmic-ray transport |[83l|8l|, while small halo models 
with L < 2 are strongly disfavored from complementary obser- 
vations I1011ll061[32l and are to be considered as extreme cases 
— we complete this discussion in Sect. 4.3.2 Therefore, we 
will stick to our median transport model which can be thought 
of as rather conservative. 

We finally confront the predicted antiproton fluxes with the 
PAMELA dat^ An error box is defined as having the en- 
ergy bin width and the 1-cr error bar half-height. A model 
is considered in tension with the data when the predicted flux 
(adding up the primary and secondary contributions) is above 
an error box without crossing i|^ A few examples are shown 
in the right panel of Fig. [7] where the same injected spectra as 
in Fig. [6] have been used an compared to the case of an anni- 
hilation into bb. A conventional annihilation cross section of 
(crv) = 3 X 10"^^ cm^/s was assumed, for which all models ap- 
pear in excess with respect to the data. This clearly illustrates 



to Eq. ( 29 ), the antiproton flux at the Earth is given by the fol- 



^ Other data are not considered because data taken from space are supposed 
to have smaller systematic errors than balloon data for a given statistics; more- 
over, PAMELA operated during low solar activity periods, which is optimal for 
low energy cosmic ray studies. 



11 




Figure 7: Left: secondary antiproton flux (background) and rough theoretical uncertainties. Right: examples of primary fluxes with the same scalar and singlino 
mass combinations as in Fig. [6] they are compared to a bb spectrum similar to the one used in | 32 1. 



the rejection power of the antiproton analysis. In the following 
we present the results we obtained from scans that cover the 
CoGeNT region, as detailed in Sect. [3] 

4.3.1. CoGeNT region scan analysis 

Fig. |4] shows the results of two scans we performed to cover 
the CoGeNT region in the plane - cr^^. Dark red points are 
those which fulfil the relic density constraint, and the light or- 
ange points are those in excess with respect to the antiproton 
data. Both populations can roughly cover the whole region, for 
small (upper panels) as well as for large values (lower panels) 
of tanyS, independently of Tqqd- Still, taking small values of 
tan 13 results in more cosmologically relevant configurations in 
average, as it can be observed. Indeed, tan/5 = 10 implies a cou- 
pling c^h a factor of ~ 10 larger than in the case of tanyS = 50, 
for a given spin-independent cross section (see Eq. [21]), which 
provides more room to get the correct relic density only from 
the annihilation into h h. Modifying the QCD phase transition 
temperature leaves these statements valid (see left/right pan- 
els). More generally, Fig.|4]demonstrates that we cannot draw a 
generic antiproton exclusion contour in this plane, since many 
configurations pass both the cosmological and the antiproton 
tests. Nevertheless, it shows that the antiproton constraints are 
relevant everywhere in this plane, which is by itself an impor- 
tant result. In the following, we provide details about the nature 
of those points which are (or not) in excess with respect to the 
antiproton data. 

In Fig. [8j we replot the scan results of Fig. [4] in the plane 
{o-v)Q-m^. Note that ((Tv)o refers to the annihilation cross sec- 
tion relevant to indirect detection, i.e. the velocity-independent 



contribution given in Eq. ( 12 ), while the full velocity-dependent 



annihilation cross section was used in the relic density calcu- 
lation. We remark that the relic density constraint obviously 
translates into an upper bound ((Tv)o ^ 5(8) x 10"^^cm^/s 
for Tqcd = 150(400) MeV, but that the cosmologically viable 
points spread down to arbitrarily small values of ((Tv)o that cor- 
respond to cases where the P-wave is dominant, i.e. where an- 
nihilation proceeds mostly into hh or a a at freeze out. This 
demonstrates that the kinematic condition 2m^>ma + mh does 
not generically ensure that annihilation into a h dominates. In- 
deed, suppressing the coupling between the singlino and the 
pseudo-scalar a is enough to suppress the S-wave, as shown 



in Eq. (12). We will come back to this later. The population 
located close to the upper cosmological bound is featured by 
those configurations for which annihilation is saturated by the 
S-wave, i.e. annihilation into ah. This is precisely where the 
antiproton constraint is the strongest. There is still a signifi- 
cant population that remains unaff'ected by the antiproton con- 
straints in this region. This actually corresponds to cases where 
the decays of a and h do not produce antiprotons efficiently 
(see discussion in Sect. |4.1|). Eventually, the impact of chang- 



^We have dismissed those data points which systematically lie below the 
secondary flux prediction by more than a standard deviation. 



ing Tqcd is a bit harder to figure out from our scan results. 
Typically, considering large values of Tqcd implies a larger 
cross section at freeze out over a broader WIMP mass range, 
i.e. when m^lxf < Tqqd (see Fig. [T]). This means that going 
from Tqcb = 150 to 400 MeV, we increase the required annihi- 
lation cross section by a factor of ~ 2 in the singlino mass range 
-4-12 GeV (see Fig.[T]). 

A way to identify the impact of the S-wave content of the 
annihilation cross section on the antiproton constraint strength 
consists in determining the S-wave to P-wave ratio, as shown 
in Fig.|9] Each panel illustrates the ratio aah/bhh versus the ratio 
(^ahl^aa for diff'erent combinations of tan (10/50 on top/bottom 
panels) and Tqcd (150/400 MeV on left/right panels). The ad- 
ditional factor of 20 applied to the ratios comes from that we 
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Figure 8: CoGeNT points in the plane (crv> (S-wave only) versus m^. Top panels: tanyS = 10. Bottom panels: tanyS = 50. Left panels: Tqcd = 150 MeV. Right 
panels: Tqcd = 400 MeV. Dark red points are cosmologically allowed, light orange points are cosmologically allowed but lead to an antiproton excess. 



want to estimate their values close to the freeze out, which oc- 
curs around Xf = m^/T ~ 20 (see Eq.fTT]). Let us denote 



SPaa = ^^Clahlbaa 
SPhh = ^^(^ah/bhh 

SPtot = 20aah/(baa-^bhh) 



(32) 



where Gah, baa, and bhh are defined in Eqs. ( |12|13|14| ), re- 
spectively. In each panel of Fig. [9] going from bottom to top 
means decreasing the relative contribution of the annihilation 
into hh, while going from left to right means decreasing the 
relative contribution of the annihilation into a a. The top left 
part of each panel, where annihilation into a a is always dom- 
inant with respect to other channels, is barely populated. In 
contrast, there is a dense wake of points going diagonally from 
the bottom right part to the bottom left part, which is charac- 



teristic of an annihilation driven by hh production; indeed, as 
long as SPhh < SPaa < 1. this channel overtops the other^ 
This implies that annihilation into a a is always subdominant in 
our sample. This is partly a consequence of requiring the Co- 
GeNT region to be fulfilled, which constrains the coupling c^h 
to be large, while c^a is left free; besides, we see that going to 
larger values of tan j3 slightly alleviates this eff'ect, as expected 



(see Eq. (22)). Another cause is less trivial: increasing Cj^a does 



not only increase annihilation into a a, but also that into ah 
(see Eq. [12]); the latter is always Icinematically allowed here, 
but not necessarily the former. 

Fig. [9] further lets us distinguish a particular area, made 



^Note that the sharp trend is reminiscent from the correlation arising be- 
tween couplings A and C^h (see Eq.[23]and discussion below); the wake would 
have been more spread without this correlation. 
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Figure 9: S-wave to P-wave ratio, aahl^aa versus aah/l>hh, as calculated from Eqs. jl2|13|14} . 



shaded in each panel, that corresponds to 5" Paa > 1 and 5" Phh > 
1. This area is the region where the S-wave saturates the full 
annihilation cross section at freeze out, and is therefore fully 
fixed by the relic density (optimal case for indirect detection). 
Naturally, the largest fraction of the models inducing an excess 
in the antiproton flux (light orange points) lies in this region. 
Nevertheless, there are still some configurations unaff'ected by 
the antiproton constraints there, which are actually those for 
which a and h are too light to decay into antiprotons (through 
quark fragmentation), or with masses >10 GeV for which the 
decay into b quarks occurs at threshold and is is then inefficient 
at producing antiprotons (-[9.5,12] GeV, see Sect. |4.1| ). In the 
latter case, usually takes rather large values because of our 
kinematical condition, which also decreases the antiproton flux 
(scaling like We note that the density of unconstrained 

configurations increases with tan JS, which comes from three ef- 



fects. First, the up/down asymmetry in the branching ratio is 
strongly dependent in tanyS, and taking a large tanyS increases 
the relative decay into r leptons compared to that into c and s 
quarks (see Eq. [6]and Fig.|5]). Second, a larger value of tanyS 
induces an increase in nth for a given spin-independent cross 
section (see Eq.|22]), and either the density of points in the mass 
range rrih -10-12 GeV increases accordingly and is pushed 
towards larger values. Third, a larger value of tanyS may also 
translate into a smaller xx h coupling, leaving more room to 
the haa coupling to set the relic density from the ^-channel 
(see Eq. 12 and the left diagram of Fig. [2]). The consequence is 
that h ^ aa opens up when masses permit, which may further 
reduce the antiproton production efficiency. 

In any case, the plots of Fig. [9] show that the antiproton con- 
straints can be generically very strong whenever S Paa ^ 1 and 
SPhh ^ 1. as expected — one can safely use Eqs. ( 12|13|14 ) to 
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Figure 10: Same as Fig.[8]after applying the cut 5Ptot > 1 (see Eq.[32). 



account for this constraint. 



In Fig. 10 which depicts the indirect detection annihilation 
cross section ((Tv)o as a function of the singHno mass m^, we 
report the same samples as in Fig. [8] after applying the cut 
SPtot > 1 (see Eq.[32]). We clearly see that this cut selects all 
points for which the relic density completely sets (crv)o, a large 
fraction of which turns out to lead to an antiproton excess. The 
remaining unconstrained cases are mostly characterized either 
by ^ 2 GeV, for which a cannot decay into antiprotons, or 
the condition < 10 GeV, mostly due to the antiproton flux 
dependence on 1 /m^ (see Eq. [3Q|). This appears more clearly 
in the plane nih-nia shown in Fig.|ll| where our results are pre- 
sented with and without the cut 5'Ptot > 1 — a rectangle is 
drawn in each panel to emphasize the most constrained mass 
range. We see that once the cut is applied, all models exhibit an 
excess with respect to the antiproton data but those with m« < 2 



GeV, and (m^ + m/^) < 20 GeV — the latter condition is equiv- 
alent to requiring m^^x < 10 GeV. 

The dependence on tanyS and Tqcd arises in all Figs.[8|[9|p^ 
and[TT] Small tan values are generically more constrained by 
the antiproton data because they lead to relatively small values 
of nih for a given spin-independent cross section (see Eq. [2T] ), 
thus favoring smaller values of which increases the antipro- 
ton flux. In this regime, the alh decay into r leptons is also 
slightly less significant than in the large tanyS regime. As to 
the impact of Tqcd, things are much simpler. The cross section 
needs to be relatively larger for < 20 Tqcd, which makes 
the indirect detection constraints more stringent for large val- 
ues of Tqcd. 

Coming back to the full results shown in Fig.[8]in the light of 
the previous discussion, we may summarize our results with a 
generic upper limit on the annihilation cross section relevant to 
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the singlino-like phenomenology as 



for 



10"^^ cmVs 



Pq 10 A GeV I cm' 



m;,/ 10 GeV 



(33) 



) Ma^mh G [2, 9.5] GeV 
UAot>0.3 



This is not an accurate Hmit, but rather a way to avoid regions 
of the parameter space that could lead to an antiproton excess. 



This is illustrated in Fig. 12 where we used the worst situation 
for the antiproton constraints, tanyS = 50 and Tqcd = 150 MeV. 
We see that a looser cut of 5'Ptot > 0-3 still leads to an excess in 
the specified mass range. 



Finally, we show in Fig. 13 how the antiproton constraint 
may translate into a constraint on the spin-independent cross 



section (see Eq. 21 ), comparing full samples (top row panels) 
to the same ones but after applying the cut 5'Ptot > 1 (mid- 
dle row panels). In each top panel, the region where the relic 
density is set by XX hh, P-wave case irrelevant to indirect 
detection, covers almost the full area but the extreme right-hand 
side, where a /z becomes dominant. The (anti)correlation 
between ((Tv)o oc [c^^a/f^x^^ o-^^ oc [c^hl{mlcp)f appears 
precisely there, which is more explicit in the middle-row pan- 
els, where the cut 5'Ptot > 1 was applied. From these plots, one 
can extract the following approximate upper limit on the direct 
detection phase space: 





2 








- K - 



< lO-^GeV" 



10-4 GeV-2 



. (34) 



In the bottom row panels of Fig. 13 we report a more ex- 
plicit translation of our antiproton constraints in terms of spin- 
independent cross section, which can be directly compared 
with Fig.|4] We have applied the cut 5'Ptot > 1- We see that 
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Figure 12: Example of generic antiproton constraints on singlino models. 
See Eq. pS). 



antiprotons provide quite strong constraints for m^^^ < 10 GeV 
in this regime, independently of tan/? and Tqcd- This nicely il- 
lustrates why the antiproton signal must be considered carefully 
in the singlino-like phenomenology. 

4.3.2. Uncertainties and potential detection prospects for 
PAMELA and AMS-02 
There are diff'erent sources of theoretical uncertainties that 
might aff'ect our predictions. Among important eff'ects which 
are related to completely diff'erent domains, some uncertainties 
may be connected to the QCD phase transition in the early uni- 



16 




Figure 13: Potential correlation between (crv>o oc [c^^^a/^x^^ andcr^^ oc [c;^./^/(m^ c^)]^. Top/bottom panels: tanyS =10/50. Left/right panels: Tqcd = 150/400 Me V. 
A cut 5Ptot > 1 was applied to all panels. 



verse (relic density), the strangeness content of nucleons (spin- 
independent cross section), the solar modulation modeling and 
the cosmic-ray transport parameters (antiproton flux). 

Uncertainties in the QCD phase transition have been shown 
to lead to a factor of 2 uncertainty in the annihilation cross sec- 
tion at freeze out in the WIMP mass range -4-12 GeV, with 
larger transition temperatures inducing an annihilation cross 
section for m^/xf < Tqqd larger by a factor of 2 than the one 
obtained for m^/xf > Tqqd, due to current constraints on the 
dark matter abundance (see Fig.[T]). We tried to account for this 
uncertainty by using both a sharp QCD phase transition and two 
difl'erent temperatures. 

Regarding uncertainties in the hadronic matrix elements, we 
refer the reader to Refs. IgOl El Ell |63l [Ml f or more details. 
The consequences on our analysis are twofold, since these un- 
certainties can lead to select either higher (larger strangeness 
fraction) or smaller (smaller strangeness fraction) values of 
Cxhlicpn?^ — see Eq. (22). In the former case, then annihi- 



much and lead to a too small dark matter abundance. On that 
account, it may become difficult to find cosmologically viable 
configurations, especially at small tanyS — this is illustrated in 



lation into h h through the t-channel would be increased, which 
would likely make the population obeying SPhh < SPaa < 1 
more prominent, thereby leading to a P-wave dominated anni- 
hilation cross section. This would further squeeze the config- 
uration range available to indirect detection. Nevertheless, we 
point out that requiring a correct relic density poses significant 
constraints then, since the annihilation rate may increase too 



the top panels of Fig. [13] where the upper regions are shown 
less favored by the relic density constraint. In contrast, the 
latter case, namely selecting smaller values of C^h/((^j3^l)r is 
more favorable to indirect detection, since a correct relic den- 
sity can be achieved by means of a larger coupling C^a^ increas- 
ing the annihilation into ah; so that the antiproton constraints 
get stronger. 

Coming back to astrophysics, the low energy cosmic-ray flux 
below a few GeV is sensitive to the solar activity. In our cal- 
culation, the solar modulation is taken into account by using 
the force-field approximation 1 1041 11Q5L which is known to 
provide a rather satisfying understanding and fit, though miss- 
ing some details, of the impact of solar activity on cosmic 
rays 1 107, JOS, i09J. We refer the reader to Ref. |109| for a 
more complete review about the current state-of-the-art. Most 
of the solar modulation models are based on numerical or ana- 
lytical solution to the Parker transport equation 1 1 10 1, the force- 
field approximation providing an analytical solution by assum- 
ing perfect spherical symmetry and a static electromagnetic 
field. More complete models lead to charge sign eff'ects depend- 
ing on the solar magnetic polarity, which reverts at each 1 1-year 
solar cycle, but they are expected to be small, at the -10-20% 
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Figure 14: Signal to secondary ratio. Top: signals are computed assuming an annihilation cross section of <crv> = 10 cm^/s. Bottom: <crv> 
solar modulation with a Fisk potential of = 600 MV. Right: = 900 MV. 
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level in terms of flux at the Earth. This roughly sizes the theoret- 
ical uncertainties afl'ecting the force-field approximation. This 
has to be compared with the error bar associated with the lowest 
energy PAMELA data points, typically ^1.5 GeV, which turn 
out to be the most constraining in our analysis (see Fig. [7]) — 
the observational error is ~ 40% in this energy range, which is 
twice more important than the solar modulation error. We have 
tuned our Fisk potential to a low/intermediate value of = 600 
MV, a rather conservative assumption with respect to the solar 
activity monitored during the PAMELA data taking fill] . This 
means that the error boxes that we have considered to derive 
our constraints safely encompass the theoretical uncertainty af- 



fecting the solar modulation modeling. In Fig. 14 we show the 
ratio of the total flux to the secondary flux, which aims at illus- 
trating the spectral break due to additional primaries for difl'er- 



ent amplitudes characterized by the annihilation cross section 

— (crv) = 10"^^cm^/s (10"^^cm^/s) on top (bottom) panels 

— and for difl'erent solar activities — = 600 MV (900 MV) 
on left (right) panels. We see that going from weak to strong 
solar activity can decrease a predicted excess down to ~ 40%. 
The PAMELA data, which were taken close to a solar mini- 
mum, are therefore very well suited for to derive low energy 
limits. 

As regards the uncertainties afl'ecting the cosmic-ray trans- 
port parameters, we stress that one of the most sensitive pa- 
rameter for dark matter-induced cosmic rays is the difl'usion 
halo size L, which characterizes the volume which the anni- 
hilation products are integrated over and which turns out to be 
strongly correlated with the difl'usion coefficient normalization 
^0 — Kq/L ^ est from stable secondary-to-primary nuclei ra- 
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tio constraints like boron-to-carbon (e.g. II112II ). Increasing L 
leads to integrate more annihilation products and thus increase 
the dark matter-induced primary flux [69, W, 113|. In this pa- 
per, we have used L = 4 kpc (see Tab.[T]). We emphasize that 
the most recent works on fitting the transport parameters tend to 
favor a larger difl'usion zone, with L ~ 5-10 kpc, which makes 
our propagation model rather conservative. These studies con- 
sider the additional input of the secondary radioactive cosmic 
ray species which allow to break the degeneracy between L 
and ^0 LSXMJ- Another unexpected problem in small-L difl'u- 
sion zone models is that they systematically lead to a secondary 
positron flux in significant excess with respect to the data at low 
energy 11141 11151 . Few other arguments in favor of larger val- 
ues of L are recalled in 1321 — see complementary arguments 
in a recent study on the diff'use Galactic radio emission 11161 . 
Our values for and L can therefore be considered as quite 
motivated and rather conservative. 

Finally, it is interesting to make some prospects for the end of 
the PAMELA campaign or the new operating space experiment 
AMS-02[^ which was successfully launched in May 2011 and 
which is supposed to take data for > 10 years. For the latter, it 
is not yet clear whether the low energy systematic errors can be 
reduced with respect to PAMELA, but statistics will for sure be 



much better. Fig.[T4]shows how a finer PAMELA/AMS-02 data 
analysis could improve our current limits or lead to the detec- 
tion of those singlino configurations which still induce cosmic 
antiprotons, but with a flux too weak to be detected or con- 
strained. This typically corresponds to an annihilation cross 
section (crv) ^ 5 x 10"^^cm^/s. We see in the bottom pan- 
els that a precision of ~ 5% makes the data quite relevant to 
such a low annihilation cross section, even in the case of strong 
solar activity. Therefore, if the PAMELA or AMS-02 collabo- 
rations manage to reach such a low precision, we expect much 
tighter constraints on light dark matter models annihilating into 
antiprotons in the future, even if the solar activity is also ex- 
pected to increase in the coming years. In particular, much 



deeper phase- space regions with 5'Ptot ^0-1 (see Eq. 32) could 
be probed. 

4.3.3. Connecting this effective approach to more constrained 
particle phenomenologies 

The question may arise of how to connect this efl'ective ap- 
proach to existing and more constrained particle phenomenolo- 
gies beyond the SM. 

First, an interesting and very general link can be made by fo- 
cusing on the kinematics of the final states. Indeed, we have 
considered that annihilation proceeds into two bodies that may 
further decay into quarks, leading to the production of antipro- 
tons. We have shown that each time the annihilation cross sec- 
tion is S-wave dominated, then this kinematical configuration 
may lead to an antiproton excess quite generically. Our eff'ec- 
tive approach is therefore fully relevant to cases for which such 
two-exotic-body production is predicted, as in e.g. Ref. I117L 
where the so-called DAMA/CoGeNT signal is explained from 



http : //www . ams02 . org/ 



a model incorporating two additional light gauge bosons that 
dark matter can annihilate into. Since one of these gauge fields 
has a mass larger than 2 GeV, it can in principle decay into an 
antiproton. This model can therefore suff'er the antiproton con- 
straints. Our constraints may also be important for completely 
diff'erent models, like in the recent study performed in 1 118]| 
of a light right-handed sneutrino dark matter that can arise in 
the frame of the NMSSM; annihilation can then indeed go into 
light pseudo-scalars without being velocity suppressed, while 
this case is not necessarily adapted to fit the CoGeNT region. 

Given our eff'ective Lagrangian of Eq. ([T]), a natural link can 
be made with singlet extensions of the MSSM, which was the 
original motivation for this work. This includes general ex- 
tensions as well as more constrained ones, like for instance 
the NMSSM, in which the dark matter candidate is a singlino- 
dominated neutralino. An example of a model of singlino dark 
matter can be found in Ref. 1231 , where the specific annihilation 
into light (pseudo)scalar Higgs bosons is considered. Their case 
(3), which corresponds to 2m^ > rUa + mh, can be directly com- 
pared with our analysis. In this specific configuration, they find 
\Cxh\ ~ \Cxa\ < \Mjkl with \Ahhh\ > \^aahV Having similar cou- 
plings between XX^^ ^i^d XX^ implies that they both have to be 
large to explain the DAMA/CoGeNT signal, which also implies 
that the S-wave contribution of the ^-channel annihilation dia- 
gram into a/z is always significant. The fact that \Ahhh\ ^ \^aah\ 
induces a competition with the ^--channel annihilation into h h 
(P-wave), but more suppressed than the ^-channel annihilation 
into a h (velocity suppression, and a squared propagator factor 
of ~ 1 /(4m^-ml)^ as compared with ~ 1 /m^ for the ^-channel). 
The antiproton limit must therefore be regarded in this case. 

For comparisons with the NMSSM, we may use Refs. 1221 
11191 11201 1121L Although the reason is not clear from the 
mass matrix elements only |^|, it appears to be difficult to 
generically find viable NMSSM configurations obeying 2m^ > 
Ma + Mh. Bosons a and h can be found very light, but not si- 
multaneously 1 11^, 11201 1121II — an exception can be found 
in Ref. [22], but the authors get an annihilation going domi- 
nantly into / /. The above cited studies still find complemen- 
tary viable regions in the NMSSM parameter space. When the 
relic density is set from the ^-channel resonant exchange of a 
light K then the annihilation is P-wave dominated and indi- 
rect detection fails to yield interesting bounds. If, in contrast, 
the main annihilation final state consists in an S-wave produc- 
tion of quarks (even if not with a branching ratio of 1), then 
an antiproton excess very likely follows in most of cases — 
this may also happen for a right-handed sneutrino dark matter 
within the NMSSM 11181 . This situation was already discussed 
in Ref. 1321 , where the antiproton limit was shown to be very 
sharp. 



5. Conclusion 

In this paper, we have designed a generic eff'ective La- 
grangian to survey the light singlino dark matter phenomenol- 
ogy, as soon as annihilation mostly proceeds into light singlet- 
dominated scalar h and/or pseudo-scalar a Higgs bosons 
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(see Sect. [2]). Within this general framework, we have com- 
puted the annihilation cross section (Sect. 2.5 ), the relic density 
(Sect. 2A\ , the spin-independent nucleon-singlino cross section 
(Sect. SfTthe antiproton spectra arising from the decays of these 
light (pseudo-)scalar bosons (Sect.|4.l|), and the subsequent an- 



tiproton flux at the Earth (Sect. 4.2.2 ). We have focused on a 
specific mass configuration that guarantees the existence of an 
S-wave contribution to the annihilation cross section, ensuring 
the relevance of indirect detection, i.e. 2m^ > nia + nih. We 
then have selected the mass and coupling parameter space to 
get a spin-independent cross section range encompassing the 
so-called CoGeNT region, and drawn several samples of ~ 10^ 
models for two diff'erent values of tanyS, 10 and 50, for which 
we have calculated the relic density. Then, we have computed 
the antiproton fluxes associated with these configurations and 



compared them to the PAMELA data (Sect. [43). 

We emphasize that there are still theoretical uncertainties af- 
fecting the estimate of the annihilation cross section at freeze 
out for ~ 10 GeV mass WIMPs, up to a factor of ~ 2 depend- 
ing on whether WIMPs decouple after or before the QCD phase 
transition. This is due to the rapid change in the relativistic de- 
grees of freedom that occurs at that time and acts as a weight 
on the annihilation cross section. We have taken this uncer- 
tainty into account by considering extreme cases, as illustrated 
in Fig.[T] 



Our results demonstrate (Sect. |4.3| ) that the antiproton con- 
straint is generically relevant to the whole CoGeNT region 
(see Fig. |4]), though there are some simple ways to escape it. 
The most trivial way to suppress the S-wave contribution to 
the annihilation cross section is to suppress the coupling be- 
tween the singlino and the light pseudo-scalar c^a^ as shown 
by Eq. ([12]). Another way is to have the light scalar and pseudo- 
scalar Higgs bosons so light that they cannot decay into an- 
tiprotons, or cannot generate antiprotons energetic enough to 
be observed. Anyway, we have shown that as soon as the 
S-wave is significant at freeze out, the antiproton signal is 
in most of the cases in excess with respect to the PAMELA 



data (see Fig. [8j which translates into Fig. [TO] after imposing 
the S-to-P wave ratio to be larger than 1). Depending on the 
singlino/scalars mass configuration, we have derived a quite 
generic upper limit on the annihilation cross section lying in the 
range {crv)^ < 10"^^cm^/s, valid for the singlino mass range 
-6-15 GeV. This limit is very competitive with respect to those 
obtained from gamma-rays |30|, and complementary to limits 
coming from high-energy neutrinos from the Sun (3T\. We re- 
call that the antiproton bounds are even more stringent for direct 
annihilation into quarks (see Fig.[7]and |[32ll ). 

Uncertainties aff'ecting our analysis have been discussed 



in Sect. 4.3.2 where we have shown that our procedure is to 
be considered as rather conservative. In addition to the points 
raised there, we remind that our antiproton flux calculation is 
poorly sensitive to the details of the dark matter distribution 
far away from the Earth, in contrast with the calculation of the 
gamma-ray flux. The main parameter which sets the flux am- 
plitude is the local density p©, so one can easily rescale our 
results and limits by a mere factor of (p^^ Ipof'. Still, we re- 
call that we have considered a dynamically constrained profile. 



and that we did not include any substructure. Although their 
survival against tidal eff'ects and encounters with baryonic sys- 
tems (stars and/or Galactic disk) still needs to be clarified, the 
survival fraction is expected to be larger and larger for smaller 
and smaller substructures, which are more concentrated II122II . 
Self-consistency would then require to include their eff'ect in 
our calculation | 97 |, which turns out to be a slight increase of 
the antiproton flux at low energy |[98l[99l. This would therefore 
make our present constraints even more stringent. 

Finally, we have shown that provided systematic errors can 
be reduced down to ~ 5%, PAMELA, from an improved anal- 
ysis, or AMS-02 could have the capability of detecting a still 
unseen signal characterized by a spectral break (see Fig. [14]) 



which would hide just within the current PAMELA error boxes 
we have considered in this paper. A null detection with bet- 
ter errors would in any case significantly ameliorate the limits 
derived in this paper. 
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Appendix A. Annihilation cross section 

In this appendix, we summarize our calculations of the an- 
nihilation cross section and of the relic density for pedagogical 
purposes (mostly dedicated to those who wish to revisit the case 
without any formal calculation software). The scattering ampli- 
tudes associated with the process ;^(/7i);^(/72) — > (pi{ki)(pjikj) 
and consistent with the Lagrangian of Eq. ([TJ, the Feynman di- 
agrams of which are shown in Fig. [2] are the following (for the 
s, t and u channels, respectively): 

Ms = ^ — ^ Kpi) C^k u{pi ) 



s - m^ -\- iTk mk 



Mu 



1 



t - m^^ 
1 



V(P2) C^j (4t + W;^) C^i U{pi) 



vipl) Cxi (4u + m^) Cxi u(pi) 



(A. la) 
(A. lb) 
(A.lc) 



where the couplings A and C are defined in Eqs. ([T]) and s, 
t and u are the Mandelstam variables and qt and qu the related 
4-momenta: 



s = (pi-\- pif = (ki + kjf 



t = q^ 



u = ql; (A.2a) 



qt = Pi- ki = kj -p2\ qu = pi- kj = ki - p2 (A.2b) 
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We also introduce the functions g(s) and h(s) such that in the 
center of mass system: 

2 



t = - g(s) + h(s) cos(0) 

u = m^^ - g(s) - h(s) cos(^) 



with g(s) = {s-m} - mj)/2 ; 



his) ^ 



s - 4m^^) 



(A.3a) 
(A.3b) 
(A.3c) 

(A.3d) 



where fij(s) is further defined in Eq. ( A.5 ), and 6 is the angle 
between pi and kf in the center of mass frame. 

The diff'erential annihilation cross section is given by: 



dcr 



JLIPS 



12^0- 



1 



sis 



2 1^1' 



■X d.o.f 



(A.4) 



4 El E2PX2 
with Pki = ^j(Pk'Pl)^ - mlm]l{EkEi) , 

where g^ is the number of degrees of freedom of particle x 
(g^ = 2 here), Sij = 1 + Sij is a symmetry factor, and where 
the integrated Lorentz invariant phase space factor reads: 



/■^'"^ = fid 



flEi {InflEj 

X{27lf6\pi+P2-ki-kj) 

Us) 



(A.5a) 



67r j_i 



with fij{s) = 



dz 



(rrii - Mjf 



(mi + rrij)^ 



1/2 



(A.5b) 



Note that z = cos(6). 

To compute the temperature-averaged annihilation cross 
section (crv), we need to express the velocity v, which 
is not the relative velocity v^/c = c\pi/Ei - ^2/^2! = 

-^jSj^ + c'^Kpi/Ei) X (^2/^2)P, but the M0ller velocity v/c = 
P\2 EH (43. Incidentally, it is useful to note that (EiE2crv) is 
Lorentz invariant. Now, let us express the cross section in terms 
of a sum of the diff'erent contributions (s-, t-, u-channels and in- 
terferences), for a deeper understanding: 

cri2/3i2 = CTuPuls-^CTuPult-^cruPulu-^ (A.6) 

O-ufiulst + Cri2^12lsu + CTuPultu • 

We are armed enough to derive the cross section in the most 
general case, i.e. we can take particles (pij^k having all a scalar 



^We stress that this general writing introduces non-physical degrees of free- 
dom when going to the total squared amplitude, notably for the interference 
terms (mixing of two different final states, e.g. scalar-scalar with pseudo- 
scalar-scalar). One has to remember that interaction terms involving c^iXc^i are 
forbidden — this is clear from the Lagrangian of Eq. Q — and can therefore be 
dropped. As a more general rule, one can replace c^iXc^j — > c^iXc^j (1 -6ij). 
Despite this caution, which is relevant when trying to give physical interpreta- 
tions to the general-case result, such a general approach proves powerful since 
a single expression can be used for diff'erent annihilation final states; one has 
just to switch the appropriate couplings on/off" to get the selected one, which 
automatically removes the unphysical terms. 



and a pseudo-scalar component; we will take care of remov- 
ing the non-physical degrees of freedom at the end^. Denoting 
crj3 = (TnPn and using g = g(s) and h = h(s), we find: 



{s-ml)^+Ylml 
K[g'-h^Y 



(A.7a) 
(A.7b) 



2h 



l-[au^2a2tg] [g^-h^]ln 



]1„ 


g-\-h 






} 




g-h 





(A.7c) 



with amt (-l)^afnt ; and index / index j for odd m 

(A.7d) 



1 KAijk(s-ml) 



4(s- ml)^ -h r; 

2o^i^? - T(o:ost + o:istg) In 
h 



g + h 



g-h 



with arnst {-^)^(^mst \ ^nd index / index j 
K 



4 



-X 



(A.7e) 
(A.7f) 



1 

gh 



{o^Qtu + o^itug ) In 



g + h 



g-h 



Parameter K, which carries dimensions of inverse squared 
mass, is defined as 



K = K(s) = 



Ms) 



6471 El E2Si 



(A.8) 



The coefficients a are associated with the numerators of terms 
involving the t- or u-channel, which all exhibit an angular de- 
pendence. This angular dependence can be expressed as a poly- 
nomial function of z reading Pxiz) = Y^n (^nxlK^) X zT. These 
coefficients are explicitly given below: 



(A.9a) 



ai, ^ -I6mlc\cl. 
a2t = -^clcj^, 



(A.9b) 
(A.9c) 



^m^i^c^i [c^j c^k(s - 4m/) -\- c^j c^k s) -(A. 10a) 
c^kclj_(g-^m^)\ 



(A. 10b) 
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4mp + c^^.m 



J) + cl;cl,m 



= 



where we have used: 



2 _ 2 



2 _ ~ ~ 



XI XJ 

(A. 11 a) 
(A. lib) 



(A. 12a) 
(A. 12b) 
(A. 12c) 



The null-velocity limit of the annihilation cross section cor- 
responds to requiring s ^^x^- ^^^^ implies that function 
h{s) 0, which a posteriori justifies the form of Eq. (|A.7|), 



from which the full limit can be readily calculated. In terms of 
the above a coefficients (we now denote them a to make ex- 
plicit that they also have to be taken in the limit s ^ A m^^), we 
obtain the following suitable expressions: 



^ ^lt+u+tu,0 
^lst+su,0 



4m2/r^2 c2 

X ijk xk 

{Aml-ml)^+Ylml 



(A. 13a) 
(A. 13b) 



2g2(4m2) 
K ^ijk [o^ost + o^osu] (4 ml - ml) 
g(4mp{(4ml-ml)^ + Tlml) ' 

(A. 13c) 



The above expressions are the ones which are relevant to indi- 
rect detection. We note that for the s-channel, only the pseudo- 
scalar exchange gives a non-zero contribution, as expected. It is 
also easy to check from these expressions that only annihilation 
into a scalar and a pseudo-scalar is non-zero in the null velocity 
limit. 

To get the Taylor expansion terms of the thermally averaged 
cross section of Eq. ( pT| ), we refer the reader to Ref. |44|. Defin- 
ing the Lorentz invariant function 



W(s) = 4Ei E2 cr/3 



e.g. from Eq. (A.4), one can demonstrate that: 



Ami 



^4mz 



(A.14) 



(A.15) 



b = 



^ml 



2W(s)-4m: 



dW(s) 



ds 



■^4mz 



The zeroth order term a (often called S-wave, though it is only 
a part of it) is therefore fully determined by Eq. ( A.13| ), while 
the first order term b requires a bit more of algebra. The lat- 
ter case is left to the patience of the reader, though explicit re- 
sults are given for the cases of XX hh and XX clcl 
Eqs. ([13]) and ([14]), respectively. For the indirect detection sig- 
nal, we have used the analytical expression obtained for the S- 
wave, while for the relic density calculation we have computed 
the thermal average numerically following Refs. L46n47il . We 



recall the full form of the thermally averaged annihilation cross 
section below: 



{CTV) = 



I lOn? rffep fx(Pl)fx(P2) Wis) 



, (A.16) 



where function / is the WIMP phase-space distribution, usu- 
ally taken as a Maxwell-Boltzmann distribution. The above 
expression can be further developed and simplified, and is fi- 
nally found to be a temperature-dependent integral over sim- 
ply two variables: s, the squared center-of-mass energy, and 
yu, the cosine of the angle between the annihilating particles — 
see ||46l|4Tl for more details. As shown above, the integral over 
[I is performed analytically in our case. In practice, we com- 
pute (crv) numerically from Eq. ( |A.16 ), where function W{s) 
is tabulated up to ^js -4m^^ = 10 m^, including all potential 
resonance and thresold eff'ects. The Taylor expansion provided 
in Eq. ( [TT) is only meant to clarify the discussion. It is indeed 
well known that such an expansion is a bad approximation in 
some cases 1 123 1. Still, the S-wave expression given in Eq. ([12]) 



turns out to be quite accurate to compute the indirect detection 
signals, which is not surprising given that the WIMP velocity 
distribution in the Galaxy peaks around x = m^/T ^ 10^ . 



Appendix B. Relic density 

We refer the reader to Refs. ||46l|43 for more details on the 
relic density calculation, which we have followed in this study. 
We still provide a few pieces of information below, mostly re- 
lated to the eff'ective degrees of freedom, which is an indepen- 
dent part, except for the last paragraph which gives some infor- 
mation about the numerical algorithm. 

Assuming entropy conservation and a Hubble rate set by the 
energy density p as ^ = 87rp/3, the Boltzmann equation that 
drives the evolution of the dark matter fluid of density n = sY 
in the early universe reads: 



dY 
dx 




Y2\ 



(B.l) 



where Mp is the Planck mass, s = (p -\- p)/3 is the entropy 
density (p being the pressure). The dynamical variable is x = 
m^/T, increasing when the universe expands. We further define 
the eff'ective relativistic degrees of freedom ggff and h^^ relevant 
to the energy density and the entropy density, respectively, such 
that: 



1/2 _ 



g.dT)T^ 



30 

^Ks(T)T' 



(B.2) 



1 -h 



T dhQff 



3/Zeff dT 
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It is useful to recall the expressions of the energy density and 
the entropy density, for a fermion or a boson species labelled /: 



Pi,hii{T) 



r cPp 



)^exp{f}Tl 
gi C (Pp E + p-v/3 



(B.3) 



J (2n 



T J (2;r)3 expjfjTl 



These integrals can actually be evaluated quite easily from se- 
ries involving modified Bessel functions of order n, Kn, which 
quickly converge to the exact results such that each species / 
carries II124II : 



?eff,/,b/f 



15 



£(±ir- 



(B.4) 



'^eff,/,b/f 



— Ki(nxi)+ -^K2(nxi) 
n 

45 
4^ 



y?. 4 x^ 

-^Ki{nxi)+ —^K2(nxi) 
n n"^ 



where Xi = rUilT. Note that in the ultra-relativistic limit Xi 0, 
these expressions tend to gi and (7/8) gi, respectively for bosons 
and fermions, as expected. We have used the following property 
of Bessel functions 11251 : 



1^ 



Y{a) Kia^ (pi y) 



Given these expressions, one can calculate the complete set 
of eff'ective degrees of freedom by summing over all physical 
degrees of freedom at a given temperature. Before the QCD 
phase transition in the very early universe down to T ~ M^, 
this amounts to all degrees of freedom of the standard model, 
including quarks and gluons, except the Z and Higgs 
bosons (the electroweak phase transition occurs earlier when 
T > Mw)- After the QCD phase transition, quarks and glu- 
ons are no longer free and remain confined into hadrons; then 
resonances have also to be included in addition to nucleons, an 
can provide extra-contribution depending on whether they are 
relativistic or not 1 126 |. 

Note that in principle, one should also add those additional 
degrees of freedom coming from non-standard particles related 
to the particular dark matter model under scrutiny. In this 
paper, the light Higgs bosons and the singlino should there- 
fore be added. Nevertheless, we can safely neglect them since 
they are also non-relativistic when singlinos decouple, around 
m^/T ~ 20 (we used rria ^ I GeV and m/^ > 4 GeV, and fo- 
cused on 5 GeV ^ < 15 GeV). Our results are illustrated 
in Fig.[T] where a sharp QCD phase transition has been consid- 
ered. Such a transition implies that gi, goes to infinity because 
of the derivative of Hq^ in Eq. ( |B.2[ ). We have suppressed this 
divergence based on that the actual transition is expected to be 
smoother 1541 . We have checked that considering this eff'ect or 
not has no significant impact on our results, at the percent level. 



As regards the numerical method, we have implemented a 
Runge-Kutta-6 scheme to track the evolution of the comov- 



ing density Y featuring Eq. (BJ_). To avoid starting the cal- 
culation too early (too small x), we first estimate the decou- 
pling temperature by implicitly solving Eq. (|B.1|), demanding 



Y(Xe) = (l+6)yeq(-^e) (with 6 < 1/2, typically) — before decou- 
pling, Y tracks its equilibrium value. We then start the calcula- 
tion at Xin < Xe, up to after the decoupling. Indeed, to achieve 
an accuracy at the percent level, we need to solve Eq. ( |B.1| ) up 
to Xend such that F(xend) > 10 5^eq(-^end), bccausc of the term 
(Y^ - Fgq). In practice, we use F(xend) = 30 Feq(-^end). Then, 
we neglect Yeq and can therefore use the solution of Eq. (|9]), 
taking x>/ = Xend and performing a Simpson integral over x in 
logarithmic steps. 
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